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Skewness as a probe of Baryon Acoustic Oscillations 



(N 

o 

(N 
O 

Q 

o 

(N 



Roman Juszkiewicz 1 ' 2 , Wojciech A. Hellwing 1 ' 3 ' 4 *, Rien van de Weygaert 5 f 

^Institute of Astronomy, University ofZielona Gora, ul. Lubuska 2, Zielona Gora, Poland 
2 Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, Warsaw, 00-716, Poland 

' ! Institute for Computational Cosmology, Department of Physics, Durham University, Science Laboratories, South Rd, Durham DH1 3LE, UK 
4 In terdisciplinary Centre for Mathematical and Computational Modeling (1CM), University of Warsaw, ul. Pawihskiego 5a, Warsaw, Poland 
5 Kapteyn Astronomical Institute, University of Groningen, P.O. Box 800, 9725LB Groningen, the Netherlands 



Accepted XXXX . Received XXXX; in original form XXXX 



o 

u 

43 

a,: 

6 
& 

(N 
> 

m 
o 

(N 



1 



ABSTRACT 

In this study we show that the skewness S3 of the cosmic density field contains a significant 
and potentially detectable and clean imprint of Baryonic Acoustic Oscillations. Although the 
BAO signal in the skewness has a lower amplitude than second order measures like the two- 
point correlation function and power spectrum, it has the advantage of a considerably lower 
sensitivity to systematic influences. Because it lacks a direct dependence on bias if this con- 
cerns simple linear bias, skewness will be considerably less beset by uncertainties due to 
galaxy bias. Also, it has a weaker sensitivity to redshift distortion effects. We use perturbation 
theory to evaluate the magnitude of the effect on the volume-average skewness, for various 
cosmological models. One important finding of our analysis is that the skewness BAO sig- 
nal occurs at smaller scales than that in second order statistics. For an LCDM spectrum with 
WMAP7 normalization, the BAO feature has a maximum wiggle amplitude of 3% and ap- 
pears at a scale of ~ 82h~ 1 Mpc. We conclude that the detection of BAO wiggles in future 
extensive galaxy surveys via the skewness of the observed galaxy distribution may provide us 
with a useful, and potentially advantageous, measure of the nature of Dark Energy. 

Key words: cosmology: theory, cosmological parameters, dark energy, large-scale structure 
of the Universe 



1 INTRODUCTION 

We live in the era of precision cosmology. The last twenty years 
of tremendous progress in astronomical observations allowed us to 
establish the standard cosmological model, the Lambda Cold Dark 
Matter model (LCDM). The cosmological parameters characteris- 
ing the model have been determined to an accuracy of better than a 
percent. 

Notwithstanding the success of these measurements, we are 
left puzzled at finding a universe whose dynamics and fate are dom- 
inated by a dark energy whose identity remains a mystery. It is not 
even sure whether it really concerns an energy component to be as- 
sociated to a new species in the Universe, a cosmological constant 
or a modification of gravity itself. It is far from trivial to constrain 
the nature of dark energy, due to the relatively weak imprint of the 
equation of state o f dark energy, in com bination with sizeable ob- 
servational errors dFrieman et al. 1 120081) . Its principal influence is 
the way in which it affects the expansion of the Universe. In turn, 
this determines the growth of the cosmic density perturbations. 

Since the discovery of the accelerated expansion of the 
universe on the basis of the luminosity distance of supernovae 
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jRiess et al] 1 19981 ; IPerlmutter et alj 1 1999b . a range of alternative 
probes have provided additional constraints on the nature of dark 
energy. The temperature perturbations in the cosmic microwave 
background, weak gravitational lensing by the large scale matter 
distribution and the structure formation growth rate are well-known 
examples of dark energy probes. 

The particular probe that we wish to address in this study are 
baryonic acoustic oscillations (BAOs), the residual leftover in the 
baryonic and dark matter mass distribution of the primordial sound 
waves in the photon-baryon plasma in the pre-recombination uni- 
verse. Following the decoupling between matter and radiation, the 
primordial sound waves devolve into slight and sub tle wiggles in 
the distribution of galaxies, the "b aryonic wiggles" l lPeebles & Yul 
ll97d : ISunvaev & Zeldovichll970l) . These have a characteristic size 
in the order of the sound horizon at decoupling or, more accu- 
rately, the baryon drag epoch at which baryons are released from 
the Compton drag of photons. The value of the latter is around 
~ 150k" 1 Mpc for the standard LCDM cosmology. 

The BAO wiggles were first detected in the power spec- 
trum of galaxies in the 2dFGRS dCole et alj l2005h and in the 
in the two-po int correlation functio n of the SDSS galaxy red- 
shift surveys l lEisenstein et alj 120051) . In th e meantime, this has 
been followed up by a flurry of studies (e.g. iTegmark et al. 20061 : 
iHutsilbOod : IPercival et 111120071: IPadmanabhan et alj|2007f) . and 
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have led to the initiation of a number of large observational pro- 
gram. Notable examples of surveys that (partially) probe dark en- 
ergy by means of Bar yonic A coustic Oscillati ons are WiggleZ 
jDrinkwater et alj2010h . BOSS dRoss et aljioich and the Dark En- 
ergy Survey, along w ith the ambitious ESA Euclid Mission project 
jLaureiis et al.l201lh . 



2 BAO IN WEAKLY NON-LINEAR PERTURBATION 
THEORY 

In this section we first outline the basic perturbation theory con- 
cepts underlying our BAO analysis, followed by a presentation of 
the formalism to incorporate the baryonic acoustic signatures in the 
power spectrum to evaluate the skewness 53 of the field. 



The basic idea behind the use of BAOs as probe of dark en- 
ergy is to use the sound horizon at decoupling as a standard ruler. 
Its physical size is known precisely from first principle, and its 
angular size was measured accurately in the angular spectrum of 
the microwave background temperature anisotropies. By assessing 
the size of the BAO characteristic scale over a range of (lower) 
redshifts, we may directly measure the angular diameter distance 
as a function of z, and hence constrain the nature of dark en- 
ergy. The BAO scale can be inferred from the imprint of the sound 
horizon in the spatial distribution of galaxies, and can be inferred 
from the correlation f unction, power spectrum or other measures 
of spatial clustering l Eisenstein et al.lTl998l ; Coorav et al. II2OO1I; 
Hu & Haimaijl2003i;lBTake & Glazebrookll2003l : ISeo & Eisensteinl 
2003l ; IZhan & KnoxEoOfj) . It was also suggested that the BAO sig- 
nature could be measured with a use of higher order statistics. How- 
ever th is mainly concentrated on bispectrum or three -point function 
alone dSefusatti etat]|20()6l ; iGaztanaga et alj|2009l) . Albeit, while 
the theoretical underpinning of BAOs is quite straightforward, the 
practical limitations and complications are considerable. A range 
of systematic effects, ranging from non-linearities to redshift space 
distortions and unknown galaxy biasing effects, severely compli- 
cates the analysis. 



To alleviate some of the practical complications, in this paper 
we advocate the idea to look for the BAO wiggles in the reduced 
third moment of the density field, the skewness S3. According to 
perturbation theory, the S3 is sensitive to second-order logarithmic 
derivative of the matter power spectrum P(k). The direct implica- 
tion of this is the imprint of the BAO wiggles in the skewness. This 
forms a potentially powerful means of probing the baryonic acous- 
tic oscillations. Unlike the more conventional second order probes, 
i.e. two-point correlation function and power spectrum, the skew- 
ness is insensitive to the bias of galaxies in case this involves simple 
linear or local bias. As important is the fact that skewness has been 
found to be rather marginally affected, in comparison to two-point 
correl ation function or power spectrum, by redshift space distor - 
tions jBouchet et all 1 9921 ; Ifflvon et al.lll995l ; lBouchet et all 19951) . 



These observations suggest that the skewness may offer a 
cleaner probe of the baryonic acoustic oscillations in the galaxy dis- 
tribution, less ridden by the systematic uncertainties that still beset 
the present BAO experiments. In this study, we present the basic 
results for measuring the signature and scale of baryonic acous- 
tic oscillations. We assess the BAO signature in the skewness of 
the cosmic matter density field, by invoking the theoretical power 
spectrum for a universe filled with a given baryonic fraction of mat- 
ter. In §[2] we present the perturbation theoretical basis of the use 
of skewness as BAO probe, including a discussion of the influence 
of non-linearities on the skewness estimates. This is followed by a 
presentation of the results in §[3] towards the feasibility of relating 
the skewness to the scale and amplitude of the BAO. Prospects and 
complications are discussed in §|4] 



2.1 Skewness and Perturbation Theory 

In our discussion on the skewness S3 of the density field, we focus 
on the volume average second and third order moments. Invoking 
the ergodic theorem|jor fair-sample hypothesis, we can express the 
volume-averaged J-point correlation function as: 



tj = Vw \ dxi...dxjW(xi)...iy(xj)£ J (xi,...,xj), (1) 

J s 

where Xi is the comoving separation vector, W(x) is a window 
function with volume 



(2) 



v w = / rfxvy(x), 



and the integral covers the entire volume S. Because of the fair- 
sample hypothesis, £j does not depend on the location x and is a 
function of the window volume Vw onlv lPeeblej J 198Cj> . 

The skewness S3 is the normalized ratio of third order to sec- 
ond order moment, defined as 



S :i (R) 



( 3 (R) _ fa 



(3) 



6 J (R) ^ 

The skewness of the cosmic density field has been assessed 
within the context of a perturbation theory analysis JPeebleslll98Cl : 
ljuszkiewicz et al. 1993; Bernardeau 1 994ah . I n the linear and quasi- 
linear regime, ljuszkiewicz et al. I ( 1 19931) and iBernardeaul dl994j) 
found that the skewness S3 of the field, smoothed by a spherical 
top-hat window filter of scale R, is given by 



S 3 (R) = j +7iCR)- 



(4) 



The factor 34/7 would have been obtained if smoothing were no 
taken into account, and was inferred by Peebles for an O = 1 cos- 
mology. The smoothing introduces a dependence on scale, through 
the logarithmic derivative 71 (R) of the variance, 

dloga 2 (R) 



71 (R) = 



(5) 



dlogR 

where n 2 (R) is the variance within the spherical top-hat window 

Wth(x), 

1 



\R) = 



2tt 2 



dk k 2 P{k)WT H {kR) 



(6) 



In terms of the effective slope n e / / of the power spectrum at scale 
R, 71 (R) may be written as 

7i( J R) = -(n e// + 3). (7) 

Several authors confirmed that equation (2) can be used as a fair 
estimator of the skewness of the top-hat filtered density field on 
linear and weakly n on-linear scales (e . g. ljuszkiewicz et alJI 19931 ; 
lBernardeau|["l994allbl: iBaugh et al.ll 19951 ; Iszapudi et alJll999h . 

1 The ergodic theorem states that the ensemble average of statistical quan- 
tities over a large number of different realizations of a given probability 
distribution is equal to the average of the same quantity over a sufficiently 
large volume, representing a statistically fair sample of the statistical pro- 
cess. 
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2.2 Skewness and Power Spectrum 

Perturbation theory may be used towards calculating the expected 
shape and amplitude of the BAO signal for a given cosmological 
model with power spectrum P(k). On the basis of a few straight- 
forward mathematical operations, one may infer that the 71 factor 
eqn. $5$ can be expressed in terms of an integral over the power 
spectrum P(k) of the density field (we present the derivation in the 
appendix lAl: 



2 J X dkP(k)j 1 (kR)kR[j {kR) - 2j 2 (kR)] 
3f~dkP(k)jl(kR) 



(8) 

where jo, j\ and the are the spherical Bessel functions of the 



To obtain insight into the beha viour of 71 as a function of cos- 
mology, we use the expressions bv lEisenstein & Hul fl998lll999l) 
to model a linear power spectrum that incorporates BAO wiggles, 



P(fe) = A ■ k n ° ■ T(kf 



(9) 



where n a is the primordial power spectrum slope, usually close 
to 1.0. The transfer function T(k) encapsulates the details of the 
cosmological model, including the baryonic wiggles, 



T(k) = f b T\k) + (1 - f b )T cdm {k) 



(10) 



where ft = f2(,/fio is the fraction that baryons contribute to the 
full matter density. The effects of the baryons are multifold. In ad- 
dition to the relatively small acoustic oscillations, it also includes 
major effects like Silk damping. To investigate the relative impor- 
tance of the baryonic oscillations, we compare the value of 53 for 
power spectra based on a baryonic transfer function that includes 
the acoustic oscillations with those of "no-wiggles" power spectra 
where these signatures are averaged out. 

The amplitude A of the power spectrum is normalized on the 
basis of the value of the (linear) as parameter at z = 0, 



k 2+ns T 2 (k)WTH(k-R T H), (11) 



for Rth = 8ft" 1 Mpc. In our study, we have adopted the LCDM 
model at z = as our reference cosmology, using the WMAP7 
best-fit values f or the cosmological parameter estimates (data 



wmap7+bao+h0, iKomatsu etalJ OOlll) ): Q h 2 = 0.134, fif,/^ = 
0.0226, Six = 0.728, as = 0.809, n 3 = 0.963, h = 0.704. 

2.3 Non-linear Influences 

The linear power spectrum expression of eqn. (|9) will provide us 
with an impression of the influence of baryonic wiggles on the 
skewness S3, via the integral expression of eqn. ((8), for any given 
set of cosmological parameters. 

However, even on scales as large as those characteristic for 
BAO wiggles, R ~ IOO/1" 1 Mpc, the effects of non-linear evolu- 
tion of the density field become noticeable. Non-linear amplitude 
growth and mode coupling will affect the BAO signal, strongly 
on non-linear scales and subtly but significantly on quasi linear 
scales. A range of studies have demonstrated that non-linearities 
do indeed represent a m ajor influenc e . On th e basis of the ana- 
lytical fitting formula by Ismith et al.l < |2003|) . iGuzik et al.l J2007h 
showed that percent-level shift of the acoustic peak location in the 
two-point correlation function ma y be expected. Sign ificant shifts 
were also found by the study of ISmith et al.l d2007l) . on the ba- 
sis of perturbation theory and on the basis of the power spectrum 
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Figure 1. BAO imprint in power spectrum. To highlight the effect of the 
acoustic oscillations, we plot the ratio of the spectrum including the BAO 
wiggles to that of the corresponding power spectrum without BAO wiggles. 
Dashed line: ratio of linear power spectra with and without BAO. Dotted 
line: ratio of non-linear halofit corrected power spectra with and without 
BAO. Dashed-dotted line: ratio of non-linear halofit corrected power spec- 
tram to lineal' "no-wiggle" spectrum. The latter highlight the effects of non- 
linearities. 



of hal os in simulations. Following a similar route, lAngulo et alj 
J2008h looked at the impact of non-linearities on the power spec- 
trum of galaxies in semi-analytical models. They conclude that per- 
ce nt level shifts cannot be exclud ed. The thorough analytical study 
by ICrocce & Scoccimarrol d2008l) has provided particularly strong 
insight into this issue, supporting the claims for a non-negligible 
influence. 

To account for the non-linear power spectrum evolution 
we have decided to use the well-known halofit formula of 
ISmithetalJd2003h . We are aware of the fact that this approach may 
not be the most accurate. However we can treat this procedure as a 
first order approximation of the non-linear effects. The full study of 
the impact of the non-linearities calls for N-body experiments and 
will be presented in the accompanying paper. 



2.4 BAO impact on Power Spectra 

Figure Q] compares the power spectra - for both linear as well as 
non-linear situations - that incorporate the BAO wiggles with those 
of "no-wiggle" ones where these are smoothed out. To compare 
these we plot the ratio of the spectrum with BAO wiggles to the 
"no-wiggle" one, so that the unity baseline represents the corre- 
sponding "no-wiggle" spectrum. 

We compare three pairs of spectra. The first pair compares 
the linear BAO wiggle spectrum to the corresponding linear "no- 
wiggle" one. For both we use the outcome from the Eisenstein & 
Hu expression eqn. JlOt . The second pair consists of the non-linear 
spectrum with BAO wiggles and the corresponding non-linear one 
with wiggles smoothed out. Both spectra are obtained from the cor- 
responding linear spectra, subsequently transformed via the halofit 
formula. The third pair compares the full non-linear, halofit cor- 
rected, spectrum with baryonic wiggles to that of the no-wiggles 
linear spectrum. 

The first pair nicely illustrates the imprint of the baryonic 
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Figure 2. Reduced skewness and non-linearities. The lines depict the ratios 
of skewness computed using spectra including BAO wiggles to that com- 
puted with 'no-wiggles' spectra. Dashed line: skewness following from lin- 
ear power spectra. Dotted line: skewness following from non-linear halofit 
corrected power spectra. 
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Figure 3. The acoustic oscillations imprint in the skewness S3 . In all cases 
we show the skewness computed using non-linear P(k) with wiggles divided 
by the skewness computed for 'no-wiggles' power spectrum. The lines de- 
pict models with different values of Q m h 2 : 0.099 - the dashed line, 0.114 
- the dotted line, 0.134 (WMAP7) - the solid line and 0.148 marked by 
dashed-dotted line. For all models the baryon density has the same value 
fl b h 2 = 0.0226. 



acoustic oscillations. The impact of non-linearities is most starkly 
exemplified by the third pair, revealing a dramatic rise of the spec- 
tral amplitude at high frequencies. However, the effect of acous- 
tic wiggles is not dramatically changed when comparing the cor- 
responding non-linear spectrum to the non-linearly evolved "no- 
wiggle" one. Its principal effects is merely a lowering of amplitude, 
of not more than a few percent and mainly visible at higher order 
harmonics. 



3 RESULTS 

Having established the basic ingredients of our study, we next turn 
to the principal results of this study. 



2.5 Skewness & Non-linearities 

On the basis of the power spectra in figure [T] we may expect that 
non-linearities will slightly affect the estimates of the skewness S3. 
The lowering of the amplitude of the BAO wiggles in the non-linear 
regime should lead to a weakening of the BAO signal in S3. 

Figure|2]plots the skewness, computed from equations (|4) and 
{SJ, for linear and non-linear power spectrum. To enhance our ap- 
preciation of the effect that the BAOs have on the result, we plot 
the ratio of the skewness obtained for a power spectrum with the 
acoustic wiggles to that of the corresponding "no-wiggle" spec- 
trum". The dashed line marks the ratio of Ss(R)/ S^ owt3S es (R) 
for the linear power spectra, the dotted line for the non-linear 
halofit-corrected spectra. 

The first observation is that we indeed find a clear imprint of 
the baryonic oscillations in the skewness S3, in both the linear and 
the non-linear situation at the 2 to 3 percent level. The amplitude of 
the BAO signal in S3 is slightly smaller for the non-linear situation 
than the linear one. Of real importance is the observation that the 
the non-linear evolution of the power spectrum does not lead to a 
shift in the scale dependence of the skewness signal. The scale at 
which we find the peak and the minimum in the BAO impact on S3 
remains the same. Also, there is no change in the scale at which the 
BAO effect changes from negative to positive. As we will discuss 
later, the latter will be of crucial significance. 



3.1 Baryonic Fraction 

Figure [3] plots the ratio of the skewness for a spectrum with and 
without wiggles, S 3 (r) /S™ wig9les (r), for a range of LCDM cos- 
mologies in which baryons constitute a different fraction of the 
non-relativistic matter content of the Universe. The baryon density 
is the same for all models, f^/i 2 = 0.0226, a value which is tightly 
cons trained by the CMB t emperature anisotropy power-spectrum 
(e.g. iKomatsu et all j201llV ). The total matter density varies from 
n m h 2 = (0.099, 0.114, 0.134, 0.148). 

In all cosmologies, the BAO skewness signal consists of a 
characteristic pattern. Starting from small scales, we observe a de- 
crease of the skewness as the scale increases, reaching a minimum 
at scales of roughly 45 — 50/i _1 Mpc. Subsequently, it climbs 
towards values larger than unity, which it reaches at a scale of 
~ 55 — GOh^ 1 Mpc, after which it increases towards a peak value 
of 2 — 4% at a scale of around 75 — SOh" 1 Mpc. Towards larger 
scales it gradually falls off towards unity. 

We immediately observe that the BAO skewness signal in- 
creases strongly as fi m is lower and baryons constitute a more 
prominent fraction of matter content of the Universe. While the am- 
plitude of the BAO skewness signal is responding systematically to 
the cosmic baryon fraction, it does not seem to exceed the ~ 3—4% 
level. In this sense, the BAO imprint is smaller that that seen in the 
corresponding power spectrum of density perturbations, where the 
wiggles may amount to a ~ 10% effect. 
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Figure 4. The dependence of the sound horizon scale (r s ) on the crossing 
scale (r c ). The boxes mark the values predicted by PT skewness values (see 
Fig-HJ. Two lines show our fit to this scaling: the linear fit (the dashed line) 
and quadratic fit (the dotted line). The values of r c where rescaled to Mpc 
using value of h = 0.704. 



Table 1. The sound horizon scale r s at the baryon drag epoch and the 
corresponding crossing scale ion S3 signal r c . All models are computed 
with fixed Qt,h 2 = 0.0226. Values of r c are obtained with use of h = 
0.704. 





r s (,z d )[Mpc] 


r c [Mpc] 


Zd 


0.099 


164.82 


87.24 


1016.04 


0.114 


159.48 


85.13 


1018.14 


0.134 


152.95 


82.24 


1020.48 


0.148 


149.1 


80.54 


1021.67 



3.2 Skewness Scale & Sound Horizon 

We also make the interesting observation that the main wiggle in 
the skewness appears at a considerably smaller scale than that of 
the BAO peak in the two-p oint correlation function fe(r) (e.g. 
ICrocce & Scoccimarrol2008l) . This goes along with the finding that 
the BAO imprint on the skewness is visible over a considerable 
larger range of scales than that in the two-point correlation func- 
tion. While for the latter, the BAO signal is noticeable on scales 
from roughly 80/i _1 Mpc to 130/i _1 Mpc, the skewness is af- 
fected over a wider ranee of 20k- 1 Mpc to lOO/i -1 Mpc. 

While the amplitude change of the skewness is closely re- 
lated to the effective baryon fraction /(, of a cosmological model, 
the scales associated with the wiggle should be correlated with 
the scale of the acoustic horizon. The scale of the acoustic hori- 
zon r 3 is a fundamental aspect of the physics of baryonic acous- 
tic oscillations. It is the physical scale of the largest acoustic 
oscillations at the epoch of recombination, and as such is the 
scale imprinted in the BAOs visible in the distribution of galax- 
ies and baryons in the Universe. Determining the apparent BAO 
scale over a wide range of redshifts, and relating it to the sound 
horizon, is the goal of a large number of curr e nt and future 
galaxy redshift surveys, {g eo & Eisensteinj 120071 ; iMeiksin et alj 
Il999t ISeo & Eisens"tdnll2003l ; IWhitel2005l) . 

To enable the exploitation of the information content of the 



skewness on the acoustic oscillations, it is therefore of fundamental 
importance to establish the connection of the sound horizon scale 
with the characteristic scales we find in the BAO feature in S3. To 
this end, we compare the sound horizon scale r 3 for a range of cos- 
mologies to the crossing scale r c at which the effect of BAO wig- 
gles on the skewness reverses from negative to positive (see fig. [2] 
and[3j. This scale is clearly a characteristic for the BAO imprint on 
skewness, and may be formally defined as the scale r c at which 

A S3 (r c ) = ^) = 1 . (12) 

o 3 (r c ) 

In other words, it is the scale at which Ag 3 (r) — 1 changes its sign 
from negative to positive. We have also checked that the r c scale 
is corresponding with the location of the inflection point of the S3 
curve near that scale. 

We have measured the characteristic crossing scale for all 
models, and related them to the correspond ing sound horizon scale s 
and redshift of the baryon drag epoch - Zd fcisenstein & Hufl999h . 
We list the measured values in table[T] To obtain a good understand- 
ing of the systematic relation between the two quantities, we have 
also evaluated their values for a few cosmological models that are 
currently disfavoured by observations, involving very low and very 
high Q. m values. In figure[4]we have plotted the values of the sound 
horizon scale r a (zd) versus the measured characteristic skewness 
crossing scale r c . The figure shows that the relation is close to lin- 
ear. 

In order to assess whether a higher order function would fit 
the relation between r c and r s better, we evaluated a first order and 
a second order fit, 

r s — a; • r c + bi (linear) , (13) 

r s = a q ■ r\ + b q ■ r c + c q (quadratic) . (14) 

For the relation shown in figure |4] we have found ai — 2.365 and 
bi = -40.92 for the linear fit, with xled = 0.99. The quadratic 
fit appears to be considerably better, and has a Xred = 0.052 at 
parameters a q = 0.0187, b q = -0.934 and c q = 95.538. 

The important conclusion from this result is that it is indeed 
possible to determine the sound horizon scale at the baryon drag 
epoch, r s (zd), if we can succeed in successfully measuring the 
characteristic skewness crossing scale r c . 



3.3 Skewness BAO signature 

As a final aspect of our analysis, in figure [5] we plot the behaviour 
of the reduced skewness as a function of (top-hat) scale R, at scales 
where the BAO feature is most prominent. The red solid line rep- 
resents the skewness, as predicted by perturbation theory, for a 
power spectrum with BAO wiggles. For comparison we also plot 
the skewness resulting from the power spectrum with a 'no- wiggle' 
transfer function. 

For the spectrum with wiggles, we find a characteristic shoul- 
der at around the crossing scale. The conclusion from this compar- 
ison is that it should indeed be possible, at the 1 — 2% level, to find 
the imprint of baryonic acoustic oscillations. In an idealised survey, 
without systematic errors and bias, it should therefore be possible 
to detect the BAO feature. 
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Figure 5. Skewness BAO feature. The red solid line represents the volume- 
averaged skewness as function of top-hat scale. It shows the characteristic 
skewness BAO feature, a mild shoulder around the characteristic crossing 
scale r c . For comparison, the dotted line represents the skewness for a "no- 
wiggle" power spectrum. The insert zooms in on the shoulder at the crossing 
scale r c . 



4 SUMMARY AND DISCUSSION 

In this report we have studied the imprint of baryonic acoustic oscil- 
lations (BAOs) on the skewness of the density field, i.e. its reduced 
third moment. The skewness of the density field has been found ear- 
lier to be a robust test of the gravitational instability mechanism on 
cosmological scale. Amon gst others, it may fun ction as a probe of 
the nature of gra vity (e.g. Hellwing et al.l 201dh. and of the initial 
conditions itself jDurrer et al.ll2000t IWhitelll999h . This has prod- 
ded us to investigate whether the skewness could also help towards 
inferring cosmological relevant information via BAOs. 

To investigate the question of the BAO sensitivity of skewness, 
we have resorted to perturbation theory. From this, we obtain an 
integral expression for the skewness, dependent on the power spec- 
trum. To model the baryonic acoust ic oscillations, we have used th e 
linear power spectra expressions of lEisenstein & Hul jl998l Tl 999h . 
To take account of the non-linearities in the evolving m atter density 
field, w e use the approximate analytical expressions of lSmith et alj 
d2003h . 

We find that the BAO skewness signal is characteristic and 
may not only be detectable in the observational reality, but may 
even offer a potentially powerful alternative to existing means of 
probing BAOs. The BAO skewness signal appears to have an am- 
plitude in the order of 3 — 4%. Interestingly, it stretches out over a 
substantial range of scales, from ~ 20 - 100/r 1 Mpc and occurs 
at much smaller scales than the BAO imprint on power spectrum 
and correlation function. Perhaps most importantly, we have estab- 
lished a strict, near-linear relation between the sound horizon scale 
at the baryon drag epoch and the scale at which the BAO skewness 
signal crosses from negative to positive. 

Even though the amplitude of the BAO skewness signal is 
somewhat lower than that of the BAO wiggles in the power spec- 
trum and correlation function, it is far less sensitive to several sys- 
tematic effects that still represent a major challenge for inferring 
cosmological parameters on the basis of BAO measurements in ob- 
servational surveys. One major advantage of the BAO skewness 
signal over that in second order measures like the two-point cor- 
relation function and power spectrum, is that it is less sensitive to 
the bias of the galaxy population with respect to the mass distribu- 



tion. If the galaxy density is a local function of the mass density, 
the relation betw een the skewness and the variance of the densit y 
field is preserved ( jjuszkiewicz et al.ll995l : lFrv &Gaztanagall993h . 
In other words, the shape of the 53 reduced skewness as a function 
of scale is preserved when one concerns linear and local biasing 
of the density field. A second, and perhaps even more prominent 
advantage, is that the skewness - in the weakly non- linear regime - 
is relatively insensitive t o redshift space dist ortions dBouchet et alj 
ll992l : lHivon et alJll995l : rBouchet et aljl995h . 

While this short publication is meant to establish the feasibil- 
ity of using skewness to probe BAOs, there are numerous issues 
and details that will be discussed and evaluated in an upcoming 
study. One major issue is the influence of non-linearities in the 
density field. Subtle non-linear effects may go beyond the scope 
of what the non-linear halofit power spectra may model. Also, the 
complexities of galaxy bias may only be adequately modelled by 
means of large N-body simulations and realistic galaxy formation 
models. On the othe r hand, the curren t observational estimates of 
g 3 from th e SDSS jRoss et al.|[2007h and 2dFGRS dCroton et alj 
I2004l2007l) are characterised by relative big errors. This suggests 
that detection of the BAO feature will be difficult. In the end, proper 
modelling of hydrodynamical and radiative processes will be es- 
sential for a truly f ull treatment of non-linear and biasing effects 
dGuillet et al]|2010h . 
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APPENDIX A: DERIVATION OF THE ANALYTICAL 
FORM FOR THE 71 (R) 

To derive expression |8]f or the 71 parameter, we proceed as follows. 
The 71 parameter is defined as 

dlogcr 2 ^) 



7i 



dlogi? 



(Al) 



with u 2 (R) = (5 2 (R)} the second moment of the density contrast 
smoothed on (comoving) scale R. For a tophat filter, the latter is 
defined as 



Sr(x) 



J 5(x> 



)W R (\x-x'\)d 3 x' , 



where the expression for the tophat filter Wr(x) is 
3/4ttR 3 , x^R, 



W R (x) 



0. 



x> R. 



(A2) 



(A3) 



Using the integral expression for the density contrast, we find the 
follo wing integral ex pression for the second moment of the density 
field JPeebledTi980l) ): 

(s%m = 

d 3 x-'d 3 a-"{5(x"*')(5(a- 7 ')>W / fl(|x - x'\)W a (\x - x""\) 

d 3 x'd 3 x"£(x' - i 7 ')W fl (aT') ■ W fl (a; 7 ') . (A4) 

It is convenient to evaluate this integral expression in Fourier space. 
To this end, we use the following Fourier conventions, 

5{x) = (2tt)- 3/2 J 8 k e lU - s d 3 k, 

5 k = (2tt)- 3/2 J 5{x)e-^ x d 3 x. (A5) 
Recasting the double integral in equation [A4l to Fourier space yields 



<r 2 {R) = (2n)- 3 J J d 3 kd 3 k'{8 k 8 k ,)e lk ' 3+lk '' s W{kR)W{k'R). / 

(A6) 

leading to the following expression in terms of the power spectrum 

P(k), 

a 2 (R) = (2tt)" 3 f d 3 kP(k)W 2 (kR) 



= 7T^ I dkk 2 P(k)W 2 (kR). 



(hi) 



where the power spectrum is defined as the Fourier transform of the 
variance, 



(6 k 6 k .) = P(k)8o(k + k').. 



(A8) 



in which Sd is the Dirac delta function. The Fourier expression for 
the tophat filter, W(kR), is given by 

W(kR) = — i- (sinfci?- k R cos kR) = -ji(x) , (A9) 
{kR)^ x 

with ji (x) the first order spherical Bessel function. 
Differentiation of cr 2 (R) by R yields 



R 



l = ±fy k P m >W(kR)^, (A10) 



dR dR 
where we have used the fact that 

w <*»> 

From this, we may immediately find the expression for 71, 
dlog a 2 {R) 



71 



d Ins; It 

2^dkP{k)k 2 {kR)W{kR)^^i 



/ °° dkP{k)k 2 W 2 {kR) 



(A 12) 
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Recalling the fact that the spherical tophat filter is directly related 
to the first order Bessel function, W(x) = -ji(x), and using the 
derivation relations for Bessel functions, we find for the derivative 
of the tophat function, 

dW 3 . , . 3 ., , . 3 . . . 1 ,. , . „ . , ... 
= -^3iW--3iW = - ? 3i(it)+- bo(x) - 2j 2 {x))} . 

(A 13) 

and inserting this relation into equation I A 121 after some manipula- 
tion immediately leads to the following expression, 

2/ °° dkPjk^jk^kR [j (kR) - 2j 2 (kR)] 
11[ ' 3f °°dkP(k)jUkR) 

(A14) 

This relation for 71 is the central expression of this paper. 
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